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P REFACE 

This  Memorandum  is  part  of  a  study  requested  by  the  Advanced 
Research  Projects  Agency  of  the  utilization,  cost  and  effectiveness 
of  employing  acoustic  sensors  as  part  of  a  U.S.,  Soviet,  or  Mainland 
China  defense  against  a  low-flying  aircraft  attack.  The  summary  re¬ 
port  of  the  study  is  given  by  the  classified  RAND  Memorandum, 
RM-5713-ARPA,  Acoustic  Surveillance  for  Defense  Against  Low-Flying 
Aircraft  (U) ,  by  R.  M.  Smith. 

The  results  presented  here  on  sound  transmission  through  the 
atmosphere  and  the  resulting  sound  intensity  at  ground  level  due  to 
an  elevated  source  should  be  of  interest  not  only  in  the  field  of 
aircraft  detection  but  also  in  the  area  of  noise  control. 
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SUMMARY 


This  Memorandum  presents  the  results  of  a  study  of  the  propaga¬ 
tion  of  sound  in  the  atmosphere  with  particular  emphasis  on  the  de¬ 
termination  of  the  sound  intensity  produced  at  ground  level  by  an 
elevated  source  such  as  a  low-flying  aircraft.  For  the  purposes  of 
this  study,  a  stratified  atmosphere  is  assumed  in  which  the  tempera¬ 
ture,  wind  velocity,  and  humidity  are  dependent  only  on  the  altitude, 

A  three-dimensional  ray-tracing  program  is  developed  which  net 
only  determines  the  ray  geometry  and  propagation  time  but  also  the 
variation  of  intensity  along  the  ray.  This  computation  of  intensity 
takes  into  account  the  focusing  effects  due  to  atmospheric  refraction 
as  well  as  the  atmospheric  attenuation  along  the  path. 

The  program  is  used  to  determine  the  ground-level  sound  intensity 
due  to  a  source  at  an  altitude  of  500  ft  radiating  100  kW  (134  HP)  of 
acoustical  power.  The  intensity  is  determined  as  a  function  of  ground 
range  and  azimuth  relative  to  the  source  for  three  representative  arc¬ 
tic  meteorological  conditions.  The  results  show  that  in  the  absence 
of  atmospheric  absorption  the  intensity  in  the  immediate  vicinity  of 
the  source  varies  inversely  as  the  square  of  the  range  from  the  source. 
However,  as  the  range  increases,  atmospheric  refraction  can  result  in 
shadow  zones  of  abnormally  low  intensity  as  well  as  regions  in  which 
focusing  occurs  with  a  consequent  enhancement  of  intensity. 

As  a  result,  it  appears  that  the  detection  range  of  an  elevated 
source  is  extremely  sensitive  to  azimuth  of  propagation  as  well  as 
source  altitude,  power,  and  frequency.  Thus  even  for  the  high  inten¬ 
sity  source  assumed  the  detection  range  might  be  as  short  as  6000  ft. 
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SYMBOLS 


a  =  sound  absorption  coefficient  (ft*) 

=  classical  sound  absorption  coefficient 

a  =  molecular  sound  absorption  coefficient 
m  r 

ar  ■  average  sound  absorption  coefficient  in  the  r^  layer 
C  c  refraction  constant 
Cy  =  specific  heat  of  air  at  constant  volume 
=  vibrational  specific  heat  of  oxygen 
D  *=  distance  measured  along  a  ray  from  the  source 
ADr  =  path  length  in  traversing  the  rth  layer 
f  =  acoustical  frequency 

H  =  ratio  of  water  molecules  to  air  molecules 

h  «  Planck's  constant 

h^  =  altitude  of  the  top  of  the  r^  layer 

I  *=  sound  intensity  with  atmospheric  attenuation 

IQ  =  sound  intensity  with  no  atmospheric  attenuation 

K  =  reciprocal  of  the  relaxation  time  of  the  first  vibrational 
state  of  oxygen 

k  ■  Boltzmann's  constant 

M  -  number  of  the  layer  containing 

M  ■  molecular  weight  of  air 
a 

■  molecular  weight  of  water 
m  *  number  of  the  layer  containing  z 

m 

n  ■  number  of  the  layer  containing  z 
n^  “  unit  vector  along  the  wave  normal 
Pq  ■  acoustical  power  of  the  source 
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differential  power  radiated  into  a  solid  angle  df! 


=  vapor  pressure  of  water  as  a  function  of  temperature 
■  radial  distance  from  the  source 


universal  gas  constant 
relative  humidity 

differential  displacement  in  the  direction  of  V“ 
cross  section  of  a  differential  flow  tube 
differential  displacement  in  the  direction  of 

c 

absolute  temperature 

propagation  time  from  source 

propagation  time  from  z  to  z 

m 

t  h 

increment  of  t(z)  in  the  r  layer 
the  smaller  of  z  and  h 

r 

the  larger  of  z  and  h  . 

m  r- 1 

still  air  sound  velocity  at  altitude  hr 
still  air  sound  velocity 

component  of  ray  velocity  in  the  xz  plane 
component  of  ray  velocity  along  the  wave  normal 
value  of  v  at  source 
expression  for  v  in  the  rC^  layer 
easterly  component  of  wind  velocity 
northerly  component  of  wind  velocity 
x  component  of  wind  velocity  at  altitude  hf 
y  component  of  wind  velocity  at  altitude  hf 
wind  velocity  vector 
x  component  of  w 
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wy  ■  y  component  of  w 
w  «■  value  of  w  at  source  altitude 

AO  X 

wyQ  *  value  of  wy  at  source  altitude 

wxr  *  expression  for  w^  in  the  rC^  layer 

wyr  -  expression  for  wy  in  the  rth  layer 

x  ■  horizontal  coordinate  of  the  ray  measured  from  the  source  in 
the  direction  of  the  initial  wave  normal  azimuth 

XA  *  mo^  fraction  of  oxygen  in  air 

x(z)  =  x  displacement  of  the  ray  from  z  to  z 

m 

Axf  »  increment  in  x(z)  in  the  rth  layer 

y  *  horizontal  coordinate  of  the  ray  measured  from  the  source  nor¬ 
mal  to  the  initial  wave  normal  azimuth 

y(z)  *  y  displacement  of  the  ray  from  z  to  z 

m 

Ayr  ■  increment  in  y(z)  in  the  rth  layer 

z  =  altitude  coordinate  of  the  ray  measured  from  ground  level 

z^  =*  maximum  altitude  along  a  ray 

zm  *  minimum  altitude  along  a  ray 

z  ■  source  altitude 
o 

ar  *  still  air  sound  velocity  gradient  in  the  rth  layer 
0  *  gradient  of  w  in  the  rth  layer 

0yr  -  gradient  of  wy  in  the  rth  layer 

T1  ■  coefficient  of  viscosity  of  air 
0q  -  wave  normal  elevation  angle  at  source 

9  «  wave  normal  elevation  angle  at  an  altitude  z 

9r  -  elevation  angle  of  the  ray  at  an  altitude  z 
9  ■  elevation  angle  of  the  wave  normal  at  an  altitude  h 

r 

U  ■  coefficient  of  heat  conductivity  of  air 
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p  =  density  of  air 
p^  =  density  of  water  vapor 
$  =  wave  normal  azimuth  angle  at  source 

=  azimuth  angle  of  the  ray  at  an  altitude  z 
dQ  *  differential  solid  angle  at  source 
<u  =  angular  velocity  of  the  wave  normal 


I. _ INTRODUCTION 


The  purpose  of  this  Memorandum  is  to  determine  the  effect  of  at¬ 
mospheric  conditions  on  the  acoustic  detection  of  low-flying  aircraft. 
There  is  no  question  that  low-flying  aircraft  can  produce  high  sound 
intensities  at  locations  on  the  ground  in  the  vicinity  of  the  flight 
path.  However,  effective  acoustic  detection  of  such  aircraft  depends 
on  how  this  ground-level  sound  intensity  drops  off  with  distance  from 
the  source.  This  Memorandum  is  primarily  concerned  with  the  effect 
of  atmospheric  conditions  on  the  Intensity  of  the  sound  reaching  the 
ground  at  a  given  position  relative  to  an  elevated  source  of  known 
acoustic  power. 

If  the  atmosphere  through  which  the  sound  is  propagated  has  a 
uniform  temperature,  no  wind,  and  no  dissipation,  the  acoustic  rays 
emanating  from  the  source  are  straight  lines,  and  the  resulting  inten¬ 
sity  at  ground  level  is  directly  proportional  to  the  acoustic  power 
of  the  source  and  inversely  proportional  to  the  square  of  the  slant 
range  measured  from  the  source.  The  existence  of  a  uniform  atmosphere 
of  the  type  postulated  above  is  highly  unlikely  and  it  seems  more  prob¬ 
able  that  the  temperature  would  decrease  with  altitude  while  wind  ve¬ 
locity  would  increase.  In  Ref.  I,  Rayleigh  shows  that  in  the  presence 
of  a  linear  gradient  of  temperature  or  wind  velocity  with  altitude,  the 
acoustic  rays  emanating  from  a  source  approximate  catenaries.  Other 
investigators  have  refined  these  results  by  assuming  an  adiabatic  lapse 
rate  of  temperature  and  a  logarithmic  increase  of  wind  velocity  with 
altitude.  However,  in  spite  of  the  curvature  of  the  rays  introduced 
by  these  gradients  of  temperature  and  wind  velocity,  the  intensity  at 
ground  level  differs  very  little  from  that  obtained  for  a  uniform  at¬ 
mosphere. 

While  such  assumptions  regarding  the  structure  of  the  atmosphere 
have  the  advantage  of  providing  closed- form  solutions  for  the  ray 
equations,  it  must  be  recognized  that  the  real  atmosphere  is  subject 
to  large  deviations  from  its  average  wind  and  temperature  structure, 
particularly  in  the  first  few  thousand  feet  of  altitude  where  terrain 
conditions  are  a  factor. 
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Since  the  present  study  is  concerned  with  sound  propagation  in 
this  lower  altitude  range,  it  is  apparent  that  the  use  of  an  average 
atmosphere  is  unsatisfactory.  Instead  it  is  assumed  that  the  atmo¬ 
sphere  is  horizontally  stratified  and  that  within  any  given  stratum 
the  temperature,  humidity,  and  both  horizontal  components  of  wind  ve¬ 
locity  have  linear  gradients  with  altitude.  In  this  way,  the  best 
available  representation  of  the  actual  atmosphere  based  on  the  most 
current  meteorological  sounding  can  be  incorporated  into  the  computa¬ 
tions  rather  than  an  idealized  average  condition. 

In  Section  II  the  refraction  law  applying  to  this  type  of  strat¬ 
ified  atmosphere  is  developed  and  the  differential  equations  for  the 
determination  of  the  ray  paths  are  formulated  in  three  dimensions. 
Given  the  initial  ray  direction  at  the  source,  these  equations  are 
solved  for  the  x,  y,  and  z  coordinates  of  the  ray  paths  as  a  function 
of  the  propagation  time.  In  addition,  a  determination  of  the  ray  den¬ 
sity  and  the  atmospheric  absorption  makes  it  possible  to  determine 
the  intensity  at  any  point  along  the  ray. 

In  Section  III  the  method  is  used  to  determine  the  ray  patterns 
and  the  ground  intensity  variations  for  several  simple  temperature 
and  wind  gradient  conditions.  Finally,  three  representative  arctic 
weather  profiles  are  selected  and  the  ground  intensity  variations  with 
range  and  azimuth  are  determined. 
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II.  ANALYSIS 


STATEMENT  OF  THE  PROBLEM 

At  a  given  observer's  position  at  ground  level  it  is  desired  to 
determine  the  sound  intensity  resulting  from  a  sound  source  of  known 
acoustical  power  located  at  a  specified  altitude,  range,  and  direction 
relative  to  the  observer.  It  is  assumed  that  the  atmosphere  through 
which  the  sound  is  propagated  is  horizontally  stratified  such  that  at 
any  point  the  temperature,  wind  velocity,  and  humidity  are  functions 
only  of  the  altitude.  In  addition  it  is  assumed  that  there  are  no 
vertical  wind  components. 

REFRACTION  LAW 

In  order  to  describe  the  nature  of  sound  propagation  in  a  medium 
of  the  sort  specified  above,  it  is  first  necessary  to  determine  the 
refraction  law  which  applies. 

In  Fig,  1,  the  unit  vector  represents  the  instantaneous  di¬ 
rection  of  the  wave  normal  at  a  given  point  of  the  sound  wave  front. 
The  vector  v  is  the  still  air  sound  velocity  at  the  above  location 
having  a  direction  along  ^  and  a  magnitude  which  depends  on  the  local 
temperature.  The  resultant  velocity  uf  the  specified  point  of  the 
wave  front  is  then  given  as  vR ,  which  is  the  vector  sum  of  v  and  the 

local  wind  velocity,  w.  Thus  the  propagation  velocity  normal  to  the 
wave  front  is  given  by 


v 

n 


v  +  (w^  •  n^) 


(1) 


Since  the  velocity  vn  varies  over  the  wave  front,  it  can  be  shown 
that  the  wave  normal  has  an  angular  velocity,  S3,  given  by 


S3 


[Vv 

n 


x  Hj] 


(2) 
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If  a  coordinate  system  is  selected  with  the  z  axis  vertical  and 
the  unit  vector  in  the  xz  plane  then 


v  *  vn. 


(3) 


w  =  w  i  +  W  1 
x  y 


n1  ■  i  cos  fi  +  fc  sin  0 


(4) 

(5) 


where  0  Is  the  inclination  of  the  wave  normal  to  the  horizontal. 

On  the  assumption  of  a  stratified  atmosphere,  the  quantities  v, 
V  and  wy  are  functions  of  z  only.  Equation  (2)  becomes 


uu 


_  /  dv  \ 

J  \  9l  +  aT" cos  e/  cos  6 


(6) 


Thus  n1  rotates  about  the  y  axis  and  always  remains  parallel  to  the 
xz  plane.  Therefore  Eq.  (6)  can  be  replaced  by  the  scalar  relation 


d0  /  dv  \ 

dt  "  ‘  \  ai  +  ar  cos  0;  cos  9 


However,  the  angular  rate  can  also  be  expressed 


as 


(7) 


d0  90  dz  90 

dt  "  97  dF  =  ^  v  stn  9 


(8) 


Combination  of  Eqs.  (7)  and  (8)  gives 


96  . 

K  v  sln  6  +  cos  0  9wv 

2  ■  +  =  o 

c os  ^  0  dz 


(9) 


which  can  be  integrated  to  give 
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- X  +  w 

COS  0  X 


(10) 


Equation  (10)  is  the  desired  refraction  law  which  specifies  the  vari¬ 
ation  of  9  along  a  given  ray  as  a  function  of  altitude.  The  value  of 
C  which  specifies  the  ray  is  determined  from  the  initial  conditions  as 


C 


v 

o 

• 

cos  xo 

o 


(11) 


where  v  ,  w  ,  and  ©  are  the  values  of  v,  w  ,  and  9  respectively  at 
o  xo  o  x 

the  source.  For  a  fixed  source  altitude,  the  dependence  of  C  on  the 
initial  wave  normal  elevation  angle,  ©  ,  is  shown  in  Fig.  2. 

The  vertical  component  of  the  ray  velocity  vR  is  given  by 


dz 

dt 


v  sin  9 


(12) 


Thus  if  ©q  is  positive  the  ray  rises  initially  and  continues  to  rise 
until  the  angle  9  reduces  to  zero.  From  Eq.  (10)  it  is  seen  that  this 
occurs  at  an  altitude  for  which  v  +  wx  Is  equal  to  C.  At  this  alti¬ 
tude  it  is  found  from  Eq.  (7)  that 

^7=-4-(v+w)<0  (13) 

dt  oz  x 

so  that  9  becomes  negative  and  from  Eq.  (12)  the  ray  moves  downward 
until  it  either  reaches  the  ground  or  reaches  an  altitude  below  the 
source  level  at  which  v  +  w^  is  again  equal  to  C.  In  this  latter  case 
it  can  be  shown  that  9  becomes  positive  and  the  ray  begins  to  rise 
again. 

This  behavior  can  be  determined  by  plotting  v  +  w^  as  a  function 

of  altitude  as  in  Fig.  3,  for  example,  where  zq  is  the  source  altitude. 

In  this  case,  if  ©  is  zero  the  value  of  C  is  equal  to  C  .  It  is  seen 
o  o 

that  the  resulting  ray  has  a  maximum  altitude  of  z^  and  a  minimum  al¬ 
titude  of  z'  and  can  never  reach  ground  level.  As  C  increases  the 
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Fig.  2 — Dependence  of  the  refraction  constant  C  on  0 

o 


Fig.  3 — Limitations  on  the  altitude  of  a  ray 
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maximum  altitude  increases  and  the  minimum  decreases  until  at  a  value 
of  the  ray  reaches  the  ground  at  grazing  incidence  and  reaches  a 
maximum  altitude  of  z^  For  values  of  C  between  and  C2>  the  max¬ 
imum  altitude  lies  between  and  z2  while  the  minimum  is  at  ground 
level.  Finally,  if  c  is  greater  than  C2,  the  maximum  altitude  is  be¬ 
yond  the  range  of  the  indicated  meteorological  data  while  the  minimum 
is  still  at  ground  level. 

RAY  DETERMINATION 

The  three  component  equations  of  the  ray  velocity  can  be  expressed 
as 


dx 

dt 


v  cos  0  +  w 

x 


(14) 


dy 

dt  *  Wy  (15) 


dz 

dt 


v  sin  0 


(16) 


The  determination  of  the  ray  involves  the  simultaneous  solution  of 
these  equations,  together  with  the  refraction  law,  Eq.  (10).  It  is 

convenient  to  select  z  as  the  independent  variable  so  that  the  result¬ 
ing  solutions  become 


x(z) 


y(*) 


V  cos  0  +  w 

_ _ X 

v  sin  0 


w 


v  sin  0 


dz 


(17) 


dz 


(18) 


I 
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t(z) 


2 

S 


dz 


v  s 


In  6 


(19) 


m 


This  form  of  the  solution  gives  the  variation  of  x,  y,  and  t  as 
the  ray  altitude,  z,  increases  from  its  minimum  value  z^  to  its  maxi¬ 
mum  value  While  these  expressions  assume  a  rising  ray,  it  will 

be  shown  that  the  x,  y,  and  t  values  measured  relative  to  the  source 
can  be  expressed  as  a  linear  combination  of  x(z),  y(z) ,  and  t(z)  eval¬ 
uated  at  z,  zq,  and  even  though  the  ray  may  pass  through  maxima  and 
minima  on  the  way. 

In  the  previous  section,  it  is  assumed  that  v,  w  ,  and  w  are 

x  y 

functions  of  z  only.  At  this  point  it  is  further  assumed  that  the 

th 

atmosphere  can  be  divided  into  horizontal  layers  such  that  in  the  r 

layer  the  quantities  v,  w  ,  and  w  have  linear  gradients  with  altitude 

x  y 

as  follows: 


v  =  V  +  or  /z  -  h  ,\ 
r  r-1  r  \  r-1/ 

w  **  W  ,  +  B  (z  -  h 
xr  x,r-l  xr  \  r-1/ 

w  -  W  +  0  (z  -  h  ) 

yr  y,r-l  yr'  r-1/ 


(20) 

(21) 

(22) 


if 


h  .  «  z  £  h 
r-1  r 


(23) 


In  the  above  equations  a^,  Bxr »  ar>d  ^yr  are  8radients  v*  wx» 

and  w  respectively  while  V  .  ,  W  .  ,  and  W  .  are  the  values  of 
y  r- i  x ,r- i  y " 1  th 

v,  w  ,  and  w  at  h  ,  the  altitude  of  the  bottom  of  the  r  layer, 
x  y  i 

With  this  specification  of  the  me teoiological  conditions,  Eqs. 
(17)  through  (19)  can  be  integrated  to  give 


« 


f 
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(30) 
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At  -  t  (U  )  -  t  (u  ) 
r  r  r  r  r 


(31) 


where  the  functions  xf(z),  yr(z),  and  tf(z)  are  defined  as  follows: 


<  (z)  =  -  —  v  tan  8  -  A  t  (z) 
r  fY  r  xr  r 


(32) 


yr(z) 


B  V 

JJLJL 


a  l^2  -  a2  ^  cos  8 

r \  xr  r  ) 


8  sin  8  +  a  (a  +8  cos  8)1  (z) 
xr  r  r  xr  r 


+  —  A  t  (z) 
a  yr  r 
r  1 


if  e2  i  a2 
xr  r 


(33) 


8  v 


yrVr  fcos  8  i  2 

.  2  [  cos  8 

3^ 

r 


](-  if 


A 

+ 

2 

a 

r 


,  8  9 

sin  -  -  cos  - 

log 

/  L 

,  8  6 

^  O'  I  (z) 
r  r 

sin  -  +  cos  - 

. 

if  B  m  ±  at 
xr  r 


(34) 


:  (z)  *  —  [log  (  1  -  SlW-9)  -  B  I  (z)l 
r  ry  L  V  cos  9  /  xr  r  J 


(35) 


where 


I  (z) 

r 


7 


,2  2  log 


8  -  a* 

xr  r 


2  2 

v  111  B  -Of  sin8-a  cos  8  -  B 

xr  r _ r _ x£ 

cos  8 


if  B2  >  O'2  (36) 

xr  r  ' 


7T 


8  in 


xr 


.  T  8  +  cr  cos  8  "| 

■1  xr  r _ 

8  cos  8  +  0- 
L  xr  r  J 


if  B2  <  o'2  (37) 

xr  r  v  ' 
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(38) 

(39) 


A 

yr 


»W  .  -  8  V  , 
r  y ,r-l  yr  r-l 


(40) 


As  Indicated  previously,  the  functions  x(z),  y(z),  and  t(z)  are 
zero  at  zm  and  assume  that  z  is  increasing.  However,  these  functions 

can  be  used  to  determine  the  x,  y,  and  t  values  measured  from  the  source 
as  follows: 

If  ®  >  0  or  z  m  z 

o  o  m 


X  -  x(z)  -  x(z  ) 

0 

if  z 

0 

5 1  s 

and  ~  >  0 

dt 

(41) 

■  2x(zm)  -  x(z) 

-  X(V 

it  p-< 

dt 

0 

(42) 

-  2x(zm)  +  x(z) 

'  X<V 

if  z  5 

m 

z  <  z  and 

0 

dz 

dF>0 

(43) 

where  z^^  -  maximum  altitude  of  the  ray. 

Thus  for  these  conditions,  Eq.  (41)  determines  the  x  displacement 
as  z  increases  from  zq  to  z^  If  ^  is  not  infinite  then  Eq.  (42)  ap¬ 
plies  as  z  decreases  from  ^  to  z^  Finally,  Eq.  (43)  is  used  as  z 

increases  from  zm  to  Zq,  thus  completing  a  full  cycle  of  the  altitude 
variation. 

If  0  <  0  or  if  z  -  z. 

o  o  M 


m 


x(zQ)  +  x(z) 


dz 

dT>  0 


s?<° 

(44) 

(45) 

if 
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*  2x(zm)  +  x(zq)  -  x(z)  if  ^  *  z  *  zq  and  —  <  0 

(46) 

In  this  case  Eq.  (44)  applies  as  z  decreases  from  z°  to  z^,  Eq. 
(45)  determines  the  value  of  x  as  z  increases  from  z^  to  and, 
finally,  if  is  not  infinite,  Eq.  (46)  is  in  effect  as  z  decreases 
from  to  zq. 

The  corresponding  expressions  for  y  and  t  are  obtained  from  Eqs . 
(41)  through  (46)  by  replacing  x  by  y  or  t. 

The  foregoing  analysis  gives  a  complete  solution  for  an  acoustic 
ray  traversing  an  atmosphere  with  a  stratified  wind  and  temperature 
structure  in  terms  of  the  height  of  the  source  and  the  azimuth  and 
elevation  angles  of  the  wave  normal  at  the  source. 

INTENSITY  DETERMINATION 

Suppose  that  the  wave  normal  corresponding  to  a  particular  ray 
has  an  elevation  angle,  fi  ,  and  an  azimuth  angle,  (zero  by  defini¬ 
tion),  at  the  source.  If  this  wave  normal  is  perturbed  by  an  amount 

d®  in  elevation  and  d$  in  azimuth,  the  resulting  perturbations  in 
o  o 

in  the  elevation  and  azimuth  angles,  0R  and  0R,  of  the  ray  velocity 
are  as  shown  in  Figs.  4a  and  4b.  In  these  figures  vq,  w^,  and  w^q 
are  the  still  air  sound  velocity  and  the  two  wind  components  at  the 
source  location.  These  perturbations,  d0R  and  d0R,  can  be  used  to 
define  a  differential  solid  angle  dfi  at  the  source  in  the  form 

dfi  *  cos  8Rd8Rd0^  (47) 

or  in  the  terms  of  ®  and  $  this  becomes 

o  o 

dQ  ■  /l  -  — —  cos  0  |  cos  ©  d©  d$  (48) 

\  v0  °  /  o  o  o 

Since  w  «  v  ,  it  is  seen  that  the  solid  angle  can  be  expressed  to 
xo  o 

a  good  approximation  as 


Fig.  4 — Perturbations  of  a  ray  in  elevation  and  azimuth 
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dO  =  cos  ©  dft  df  (49) 

ooo 

If  the  source  Is  omnidirectional  and  has  an  acoustic  power  of  Pq, 
then  the  energy  flow  rate  from  the  solid  angle  dQ  is  given  by 

P 

dP  =  cos  ©  d©  d$  (50) 

4tt  ooo 

All  of  the  energy  flow  from  dfl  remains  confined  to  a  tube  bounded  by 
the  family  of  acoustic  rays  which  initially  bounded  the  solid  angle 
dQ.  Figures  4a  and  4b  show  the  vertical  and  horizontal  projections 
of  this  tube.  Since  no  energy  passes  through  the  sides  of  the  tube 
and  it  is  assumed  that  there  is  no  attenuation,  the  energy  flow  rate 
through  any  cross  section  must  be  equal  to  dP.  Sound  intensity  is 
defined  as  the  energy  flow  per  second  per  unit  of  area  normal  to  the 
direction  of  propagation,  thus 


dP  *  I  dS 
o 


(51) 


where  I  is  the  sound  intensity  and  dS  is  the  cross-sectional  area  of 
the  flow  tube.  Combination  of  Eqs.  (50)  and  (51)  gives 


P  cos  ©  d©  d$ 
o  ooo 

o  4tt  dS 


(52) 


The  cross  section  dS  is  the  product  of  dr,  the  width  of  the  tube  in 
the  vertical  plane,  and  ds ,  the  width  in  the  horizontal  plane.  Since 
dr  depends  only  on  d©^  and  ds  on  d$^,  Eq.  (52)  can  be  rewritten  as 


P  cos  © 


(53) 
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Thus  it  is  seen  that  the  sound  intensity  varies  inversely  as  the  area 
of  the  flow  tube.  For  a  constant  temperature  atmosphere  with  no  wind 
the  acoustic  rays  are  straight  lines  radiating  from  the  source  and 


(54) 


dT  -  R  cos  ®0  (55) 

o 

where  R  is  the  radial  distance  from  the  source.  Equation  (53)  reduces 
to 


I 

o 


P 

_  o 

2 

4n 


(56) 


the  familiar  inverse  square  variation  of  intensity  with  distance. 

The  equivalents  to  Eqs.  (54)  and  (55)  when  wind  and  temperature 
gradients  are  present  can  be  obtained  as  follows.  Equations  (41) 
through  (46)  define  a  family  of  acoustic  rays  in  the  xz  plane  with 
0q  constant  for  a  given  ray.  Thus 


dr  _ 1_ 

d©  “  Tv©“ 
0  0 


(57) 


where 


70 

o 


(58) 


An  examination  of  Eqs.  (41)  through  (46)  shows  that  the  general 
form  of  the  ray  equations  is 


I 


-18- 


*  ■  *  *<z)  + 


(59) 


where  the  function  F  is  a  linear  combination  of  x(z  )  x(z  ),  and 

o  m 

x(z^).  From  Eqs.  (17)  and  (59)  it  is  seen  that 


9x  5x(z) 

dz  Sz 


± 


v  cos  9  +  w 

X 


v  sin  9 


(60) 


and 


dx 

§0 


± 


9x(z) 

90 


+ 


9f(z  ,z 
o  m 


,ZM) 


a© 

o 


(61) 


Thus  Eq.  (57)  can  be  written  as 


dr_ 

d® 


v  sin  9 


± 


Sx(z)  SF 

aw  +  a© 

o  o 


(62) 


where  is  the  component  of  the  ray  velocity  in  the  xz  plane  and  is 
given  by 


V 


2  2 

v  +  w  +  2v  w  cos  9 

X  X 


(63) 


The  derivative  of  x(z)  with  respect  to  0q  can  be  expressed  by 
means  of  Eq.  (11)  in  the  form 


a© 


o 


v  sin  0 
o _ o 

2 

cos  0 
o 


X(z) 


(64) 


where  from  Eq.  (24) 
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and 


X(z) 


r*m 


(65) 


dAx 
_ r 

SC 


f(C  -  W  )  o'  +  V  0  1  "  sin  0(U  )  if 

L  xm  m  mxmJL  m  J 


r  *  m 


[<c  -  V>  “r  +  Vxr]  [sin  0(ur)  •  sm  0(Ur)] 

if  r  ^  m  and  (C-W  )  <y  +  VB  1^0 

L  xr  r  r  xrj 


c> 

r 

r  1  i  "| 

(„2  -  B2  ) 

3/2 

_v(Ur)  v(ur)J 

If  [<C  -  Wxr)  «  +  V^J  .  0 


if  r  «  M  and 


(66) 


(67) 


(68) 


(69) 


Since  the  function  F  is  a  linear  combination  of  x(z  ),  x(z  ),  and 

o  m 

X^ZM  *  ^ts  Part^al  derivative  in  Eq.  (62)  is  a  similar  linear  combina¬ 
tion  of  X(z^) ,  X(z^)  ,  and  X(z^).  These  three  quantities  can  be  deter¬ 
mined  from  the  definition  of  X(z)  above.  Thus  Eq.  (62)  can  be  ex¬ 
pressed  as  follows: 


I 


If  ©  >  0  or  if  z  ■  z 

o  o  m 


,  vv  sin  9  sin  © 

^ - 2 - ^  [X<z>  -  X<Zo>] 

o  vn  cos  0 

R  o 


if  z  ^  z  £ 
o 


*m  *"d  17  >  0 


w  sin  9  sin  © 
o _ I 

2 

vn  cos  © 

K  o 


2  P«*h>  '  x(2)  -  x<z0>]  lf37<0 


w  sin  9  sin  © 

-2 - 5 - 2  [2X<V  +  x<z>  -  x<z0>] 

v_  cos  0 
R  o 


if  z  s  z  <  z  and  *r~  >  0 
m  o  dt 


If  ©  <  0  or  2 
o 


o"  ZM 


.  w  sin  9  sin  © 

or  o _ o 

d©  2  n 

o  vn  cos  0 

R  o 


[X(2  )  -  X(2)] 


if  2  i  2  >  z  and  —  <  0 
o  m  dt 


vv  sin  9  sin  0 

—2 - j ^  CX(zo)  +  X(z)l  lf  if  >  0  <74> 

vn  cos  © 

R  o 


w  sin  9  sin  0 

—2 - 5 - °  [2X(zM)  +  X(zo)  •  X(z)] 

v_  cos  0 
R  o 


If  2,,  £  z  2  z  and  —  <  0 
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It  should  be  noted  that  Eqs .  (70)  through  (75)  become  indetermi¬ 
nate  if  either  0  or  ®q  is  zero.  However,  this  can  be  resolved  as 
follows.  If  z  approaches  z,^  in  either  Eq.  (70)  or  (71),  the  brack-  ' 

eted  terms  are  dominated  by  the  Mth  term  in  X(z)  ,  which  from  Eq.  (67) 
is 

^ .  c  r _ i _ M 

90  <c  •  +  VxM  LSin  Vl  '  sln  ®J  <76> 

and  as  9  becomes  small 

_ C 

90  =  "  V^QM  +  BxM^  sin  0 

Thus  either  Eq.  (70)  or  (71)  approaches  a  limit  of 


v  sin© 
o _ o 


%  +  0xM)  Cos2 


© 

o 


(78) 


as  z  approaches  z^.  In  a  similar  manner  it  can  be  shown  that  as  z  ap 
proaches  z^  Eqs.  (71)  and  (72)  approach  a  limit  of 


v  sin 
o 


(»  + 

m 


B  ) 

xm 


o 

o _ 

2  „ 
cos  0 

o 


(79) 


It  can  also  be  shown  that  Eqs.  (73)  and  (74)  approach  the  limit 
given  by  Eq.  (70)  as  z  approaches  z^ ,  while  Eqs.  (74)  and  (75)  approach 
the  limit  given  in  Eq.  (78)  as  s  approaches  z^ 

If  ZQ  approaches  Zm  and  0q  approaches  zero,  Eq.  (70)  reduces  to 


dr 

d© 


Cv  s in  0 


VR 


+  8  ) 

xm 


(80) 


as  does  Eq.  (71). 
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Similarly  as  z^  approaches  and  9^  approaches  zero, 
(73)  and  (74)  reduce  to 


both  tqs , 


dr 

d 


Cv  sin  9 
vR(5M  +  V 


(81) 


The  other  derivative  in  the  denominator  of  Eq.  (53)  can  be  expressed 
as 


Since 


ds  _  1 

d$  =  | 

o  o 


and 


(82) 


(83) 


(84) 


(85) 


where 
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and 


(86) 


(87) 


Since  x  is  a  linear  combination  of  x(z)  ,  x(zq)  ,  and  x^)  ,  and  y  is  a 
linear  combination  of  y(z),  y(zQ)  ,  and  y(^)  ,  their  partial  derivatives 
in  Eq.  (85)  are  similar  linear  combinations  of  the  partials  of  x(z)  and 

y(2)  eValuated  at  z.  V  and  ^  If  Eqs .  (17)  and  (18)  for  x(z)  and 
y(z)  are  generalized  to  include  f  as  follows: 


then 


x(z) 


y(z) 


v  cos  9  cos  f  +  w 

_ O  X 

v  s  in  9 


dz 


v  cos  9  s in  $  +  wy 

o _ J 

v  s  in  9 


dz 


cos  9  sin  $ 


(88) 


(89) 


sin  9 


dz 


(90) 
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and 


M  -S' 


cos  B  cos  $ 

_ c 

sin  9 


dz 


(91) 


z  z 


m 


When  is  set  equal  to  z^ro,  Eqs.  (90)  and  (91)  become 


(t^i  -o 

\  O 


(92) 


/ ^y(z) \  f  cos  6  dz 

l  )  •*«  9 

z  z 

m 


Combination  of  Eqs.  (88)  and  (93)  gives 


(93) 


(94) 


Since  «  v,  the  integral  term  can  be  neglected  and  if  Eq.  (91)  is 
substituted  into  the  expressions  for  y  corresponding  to  Eqs.  (41) 
through  (46) ,  it  is  found  that 


9$ 


x 


Thus  Eq.  (53)  for  the  sound  intensity  becomes 


P  cos  0 


(95) 


(96) 
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ATTENUATION 

Dissipation  of  energy  in  the  medium  results  in  an  exponential  de¬ 
crease  of  intensity  with  path  length  of  the  form 


I 


t  ‘aD 
I  e 

o 


(97) 


where  Iq  is  the  intensity  at  a  given  location  in  the  absence  of  any 
energy  dissipation,  D  is  the  path  length  in  tho  medium,  and  a  is  the 
sound  absorption  coefficient. 

The  principal  causes  of  sound  absorption  in  the  atmosphere  are 
viscosity,  heat  conductivity,  and  molecular  absorption.  Of  these,  the 
first  two  are  ordinarily  called  classical  absorption,  for  which  the  ab¬ 
sorption  coefficient,  a^ ,  can  be  expressed  as 


a 

c 


2  2  r 
4rrr  471 

- T~  i  + 

pv  L 


pM  R 

_ w  o 

cv(cv  +  V. 


(98) 


where 

f  -  frequency  of  the  sound  wave  (cycles/sec) 
p  -  atmospheric  density  (gra/cm3) 
v  -  sound  velocity  (cm/sec) 

"H  ■  coefficient  of  viscosity  of  air  (poise) 

U  -  coefficient  of  heat  conductivity  of  air  (cal/gm/°C) 

^  -  molecular  weight  of  air  (gm/mol) 

Rq  -  gas  constant  (cal /mol) 

Cy  "  specific  heat  of  air  at  constant  volume  (cal/mol/°C) 
ac  “  absorption  coefficient  (l/cm) 

Molecular  absorption  is  due  to  the  excitation  by  the  sound  wave 
of  the  vibrational  states  of  the  oxygen  molecules  in  the  atmosphere. 
The  absorption  coefficient  in  this  case  is  given  by 
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a 

m 


2nR  x.C . 
o  A  i 


Kf 


C  C.,(C„  +  R  )  ff2  2 
o  V  V  o  (f  +  K  ) 


(99) 


where 

x  =  mol  fraction  of  oxygen  in  air 

A 

=  vibrational  specific  heat  of  oxygen  (cal/mol/°C) 

K  =  reciprocal  of  the  relaxation  time  of  the  first  vibrational 
state  of  oxygen  (1/sec) 

a  =  molecular  absorption  coefficient  (1/cm) 
m 

The  vibrational  specific  heat  can  be  determined  as  a  function  of 
temperature  from  the  expression  (Eq.  (263)  of  Ref.  2). 


(100) 


where  T  is  the  absolute  temperature  and  from  spectroscopic  data  hv/k 
is  equal  to  2228  (Table  2  of  Ref.  2). 

In  Refs.  3  and  A,  it  is  shown  that  the  value  of  the  relaxation 
time  of  the  first  vibrational  mode  of  oxygen  depends  on  the  amount  of 
water  vapor  present.  As  a  result  the  parameter  K  is  given  by  the  re¬ 
lation 


K  -  5.758  X  108H2 


(101) 


where  H  is  the  ratio  of  the  number  of  water  molecules  to  the  number 
of  air  molecules  and  is  given  by 


M  p  T 
a*w 

M  p  T 
w*o  o 


H 


(102) 
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wherc 

M  =  molecular  weight  of  air 

a 

Mw  *  mo^ecular  weight  of  water 

Pw  *  density  of  water  vapor 

Pq  =  air  density  at  0°C 

T  =  27  3°K 
o 

T  *  absolute  temperature 
Finally,  can  be  expressed  as 

Vhp(T) 

Pw  R  T 


(103) 


where 


r^  =  relative  humidity 

p(T)  =  water  vapor  pressure  as  a  function  of  temperature 

Thus  the  absorption  coefficient  or  in  Eq.  (97)  is  the  sum  of  a 
and  arm.  However,  for  the  stratified  atmosphere  assumed  previously, 
a  varies  with  altitude  because  of  its  dependence  on  both  temperature 
and  humidity;  thus  a  more  accurate  representation  of  Eq.  (97)  is 

D 

-I  odD 

T  T  -  -A 

1  =  Ioe  "  Ioe  (104) 

The  integral  in  the  exponent  can  be  approximated  if  it  is  as¬ 
sumed  that  in  the  r  layer  of  the  atmosphere  the  absorption  coeffi¬ 
cient  has  a  constant  value,  af,  which  corresponds  to  the  average  tem¬ 
perature  and  humidity  for  that  layer.  In  addition,  it  is  assumed 

that  the  length  of  the  ray  path  in  the  rth  layer  can  be  approximated 
as 
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AD  ■=  J Ax2  +  (U  -  u  )‘ 
r  f  r  r  r 


(105) 


If  a  function  A(z)  is  defined  as 


n 

A(z)  *=  N  ’  a  AD 
/  j  r  r 

r*=m 


(106) 


Then  the  integral  in  Eq.  (104)  can  be  expressed  as  follows. 
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If  Eqs .  (96)  and  (104)  are  combined,  the  resulting  expression  for 
the  sound  intensity  becomes 
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COMPUTER  PROGRAM 

The  foregoing  analysis  has  been  programmed  for  the  JOSS*  computer, 
a  remote  console,  time-sharing  computing  system  developed  at  The  RAND 
Corporation,  In  this  program  the  stratified  atmosphere  is  specified  by 
inputting  the  altitudes  of  the  boundaries  of  the  strata  and  the  values 
of  the  still  air  sound  velocity,  the  north  and  east  components  of  wind 
velocity,  and  the  relative  humidity  at  each  of  these  altitudes.  The 
initial  conditions  for  a  given  ray  include  the  specification  of  the 
source  altitude,  the  initial  azimuth  and  elevation  angles  of  the  wave 
normal,  and  the  acoustic  power  of  the  source. 

Given  the  above  inputs,  the  computer  determines  the  x  and  y  co¬ 
ordinates  of  the  ray  and  the  propagation  time,  t,  for  specified  equal 
increments  in  ray  altitude.  In  addition,  at  each  point  the  sound  in¬ 
tensity  in  dB  relative  to  10  16  watts/cm2  is  computed  for  frequencies 
of  125,  250,  500,  and  1000  cps.  This  procedure  continues  until  the 
ray  either  reaches  the  ground  or  the  maximum  altitude  for  which  mete¬ 
orological  data  are  supplied. 

While  the  program  In  its  general  form  provides  a  complete  three- 
dimensional  solution  for  the  geometry  of  the  ray  field  as  well  as  the 
intensity  variations,  it  was  soon  found  that  for  the  purposes  of  this 
Memorandum  the  program  could  be  simplified  by  neglecting  the  lateral 
or  y  displacement  of  the  ray  since  it  has  very  little  effect  on  the 
variation  of  intensity  with  distance  at  ground  level.  Thus  the  pro¬ 
gram  was  split  into  two  parts,  one  of  which  determines  the  x  coordi¬ 
nates  of  the  ray  as  well  as  the  propagation  time  as  a  function  of  the 
ray  altitude,  z,  while  the  second  part  computes  the  sound  intensity 
at  ground  level  as  a  function  of  horizontal  distance  from  the  source. 
This  simplification  shortens  the  solution  time  considerably  with  a 
negligible  effect  on  the  accuracy  of  the  results. 
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III.  RESULTS 

It  is  now  of  interest  to  apply  the  foregoing  analysis  under  var¬ 
ious  meteorological  conditions.  For  all  of  the  cases  considered  in 

* 

this  section  a  source  with  an  acoustical  power  of  100  kW  (134  HP)  at 
an  altitude  of  500  ft  is  assumed. 

BASIC  EXAMPLES 

Before  considering  any  actual  meteorological  profiles,  it  is  use¬ 
ful  to  examine  a  few  simple  cases. 

Case  1 

In  this  case  the  meteorological  conditions  assumed  are  shown 
below. 
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W„  and  W_  are  the  northerly  and  easterly  wind  components.  Under  this 
condition  of  a  negative  gradient  of  sound  velocity  and  no  wind,  the 
rays  emanating  from  the  assumed  source  are  as  shown  in  Fig.  5,  where 
it  is  seen  that  all  of  the  rays  are  concave  upward  due  to  the  negative 
gradient  of  sound  velocity.  In  particular,  for  an  initial  wave  normal 
elevation  of  -4.73  deg  the  corresponding  sound  ray  is  tangent  to  the 
ground  at  a  range  12,100  ft  and  no  rays  reach  the  ground  beyond  this 
range,  thus  creating  a  shadow  zone.  In  Fig.  6  the  corresponding  vari¬ 
ation  of  sound  intensity  with  range  is  shown.  It  is  seen  that  the 
drop-off  is  still  essentially  an  inverse-square- law  variation  with 
range  out  to  12,100  ft.  Beyond  this  point,  the  ray  theory  indicates 


This  acoustic  power  level  is  representative  of  a  B-52  with  a 
thrust  level  of  10,000  lb  and  a  mass  flow  rate  of  180  lb/sec  for  each 
of  its  8  engines  (see  Fig.  33.42  of  Ref.  5). 
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zero  intensity;  however,  if  the  effects  of  diffraction  and  scattering 
are  taken  into  account  (see  Refs.  6  and  7)  the  intensity  in  the  shadow 
region,  while  no  longer  zero,  still  drops  off  much  faster  than  an  in- 
verse-square- law  variation. 

It  should  be  noted  that  the  patterns  shown  in  Figs.  5  and  6  are 
independent  of  the  azimuth  of  propagation. 

Case  2 

In  this  case  the  meteorological  conditions  assumed  are  shown 
below. 
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Under  these  conditions  of  a  positive  sound  velocity  gradient  with 
altitude  and  again  no  wind,  the  resulting  ray  pattern  is  as  shown  in 
Fig.  7.  Due  to  the  positive  gradient  of  sound  velocity  all  of  the 
rays  are  concave  toward  the  earth.  Thus  no  shadow  zone  is  produced 
and  long-range  propagation  should  be  expected.  The  variation  of  in¬ 
tensity  with  distance  for  this  case  is  shown  in  Fig.  8  and  again  dif¬ 
fers  very  little  from  the  usual  inverse-square-law  decrease.  As  in 
the  previous  example,  both  Figs.  7  and  8  are  independent  of  the  azi¬ 
muth  of  propagation. 

In  this  case  it  is  entirely  possible  that  sound  energy  might  reach 
a  give.,  observer  along  a  number  of  different  ray  paths  involving  one  or 
more  reflections  off  the  earth's  surface.  Under  this  condition  of  mul¬ 
tipath  transmission,  known  as  surface  ducting,  it  can  be  shown  that  at 
long  ranges  the  intensity  varies  inversely  as  the  first  power  of  the 
range  instead  of  the  usual  inverse-square  law.  However,  for  the  ele¬ 
vated  source  specified,  the  first  multiple  arrival  would  be  at  a  range 
of  about  36,000  ft  and  it  is  probable  that  the  inverse-first-power 
variation  cf  intensity  would  not  appear  short  of  100,000  ft  range. 

Thus  it  appears  that  for  the  ranges  considered  here  ducting  has  rela¬ 
tively  little  effect  on  the  intensity,  as  shown  in  Fig.  8. 
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Case  3 


The  meteorological  conditions  in  this  case  are  shown  below. 
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It  is  seen  that  the  gradient  of  sound  velocity  with  altitude  is  nega¬ 
tive,  while  the  gradient  of  W£  with  altitude  is  positive.  With  the 
introduction  of  wind  the  ray  and  intensity  variations  are  no  longer 
independent  of  azimuth.  In  Fig.  9  plots  are  shown  of  the  sum  of  the 
sound  velocity  v  and  w  for  upwind  and  downwind  propagation.  It  is 
seen  that  the  gradient  of  v  +  w^  is  positive  for  downwind  propagation 
and  negative  in  the  upwind  case.  Thus  the  resulting  ray  pattern  shown 
in  Fig.  10  is  concave  toward  the  earth  in  the  right  half  plane  and  con¬ 
cave  upward  in  the  left  half  plane.  In  Fig.  11  the  variat.ci  in  inten¬ 
sity  with  distance  for  easterly  or  downwind  propagation  is  shown  to 
differ  very  little  from  an  inverse-square  variation.  This  is  also  true 
in  Fig.  12  out  to  a  distance  of  6963  ft,  beyond  which  a  shadow  zone  is 

formed,  similar  to  that  of  Fig.  5. 

In  the  downwind  direction  for  which  v  +  wx  has  a  positive  gradient 
there  is  the  possibility  of  surface  ducting  similar  to  that  described 
in  Case  2.  However,  as  before,  any  modification  of  the  intensity  vari¬ 
ation  would  be  at  extremely  long  range. 

Case  4 

In  the  cases  considered  thus  far,  a  single-layer  atmosphere  has 
been  assumed  with  uniform  gradients  of  both  sound  velocity  and  wind 
velocity  with  altitude.  In  the  present  case  a  two-layer  atmosphere  is 
postulated  as  specified  below. 
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If  the  Initial  wave  normal  is  in  an  eastward  direction,  then  the 
x  axis  is  also  to  the  east  and  W)<  is  equal  to  W£.  Figure  13  shows 
that  v  +  wx  has  a  gradient  of  17.5  ft/sec  per  1000  ft  between  ground 
level  and  1000  ft,  while  from  1000  to  2000  ft,  this  gradient  increases 
to  37.5  ft/sec  per  1000  ft. 

The  resulting  ray  pattern  shown  in  Fig.  14  for  eastward  propaga¬ 
tion  again  shows  that  for  the  positive  gradient  of  v  +  w  the  rays 
are  concave  toward  the  earth.  However,  as  the  initial  wave  normal 
increases  from  -25  deg  to  +7.21  deg  the  ground  range  of  the  corre¬ 
sponding  rays  increases  from  1057  to  19,361  ft. 


At  7.21  deg  elevation  the  resulting  ray  is  tangent  to  the  bound¬ 
ary  between  the  two  atmospheric  layers  assumed  and  for  elevation  an- 
g  es  greater  than  7.21  deg  the  ray  spends  part  of  its  time  in  the 
upper  layer  where  its  curvature  is  greater  due  to  the  higher  gradient 
of  v  +  wx.  Thus  for  angles  between  7.21  and  8.95  deg  the  ground  range 
actually  decreases  from  19,361  to  16,458  ft.  Finally,  for  elevation 
angles  from  8.95  up  to  16.34  deg  the  range  again  increases  from  16,458 
to  21,184  ft.  At  the  latter  elevation  angle  the  ray  reaches  the  maxi¬ 
mum  altitude  for  which  meteorological  data  are  given,  and  rays  corre¬ 
sponding  to  higher  elevation  angles  are  terminated  at  the  2000-ft 
1  e  ve  1 . 

An  examination  of  the  rays  with  elevation  angles  between  7.21 
and  8.95  deg  where  ground  range  is  decreasing  shows  that  they  form  an 
envelope  or  caustic  curve  which  starts  at  the  maximum  of  the  7.21-deg 
curve  (1000  ft  altitude  and  8036  ft  downrange)  and  terminates  at 
ground  level  16,458  ft  downrange.  At  any  point  on  this  caustic  curve 
the  value  of  the  derivative  dr/dfiQ  in  Eq.  (113)  vanishes  so  that  this 
line  becomes  a  locus  of  infinite  intensity.  Actually  such  a  conclu¬ 
sion  goes  beyond  the  validity  of  ray  theory  and  in  order  to  determine 
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— Ray  pattern  (Case  4) 


the  intensity  close  to  such  a  caustic  curve  it  would  be  necessary  to 
consider  the  diffraction  effects  associated  with  a  wave  solution  to 
this  problem.  In  Fig.  15  the  intensity  variation  with  distance  based 
on  the  ray  theory  is  shown  for  the  above  case.  The  arrows  along  the 
curves  indicate  the  direction  of  increasing  initial  elevation  angle 
while  critical  values  of  the  elevation  angle  are  indicated  along  the 
curve.  Actually  the  curve  has  been  extended  to  an  elevation  angle  of 
40.5  deg  by  assuming  the  same  gradient  above  2000  ft  as  exists  between 
1000  and  2000  ft. 

It  is  seen  that  for  ground  ranges  out  to  16,458  ft  there  is  a 
single  ray  reaching  each  point  and  the  intensity  has  an  inverse-square- 
law  variation.  For  ranges  greater  than  19,361  ft  there  is  also  only  a 
single  arrival  but  the  intensity  decreases  at  a  rate  faster  than  the 
inverse  square  and  drops  below  the  curve  for  a  homogeneous  atmosphere 
at  a  range  of  about  26,000  ft.  In  the  range  interval  between  16,458 
and  19,361  ft  there  are  three  rays  which  reach  any  point  on  the  ground 
corresponding  to  the  three  branches  of  the  intensity  curve.  If  it  is 
assumed  that  the  pressure  variations  in  the  sound  wave  are  random, 
then  the  resultant  intensity  can  be  expressed  as 
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(114) 


where  1^,  I2.  and  I^  are  the  intensities  in  dB  corresponding  to  the 
three  branches  of  the  intensity  curve  for  a  given  range  while  I  is 
the  resultant  Intensity  in  dB.  Figure  16  is  an  enlargement  of  the 
multiple  arrival  region  of  Fig.  15  showing  the  three  component  curves 
and  the  resultant  intensity  curve  computed  from  Eq.  (114).  The  cusp 
in  this  curve  occurs  where  the  caustic  curve  discussed  earlier  inter¬ 
sects  the  ground. 

As  indicated  previously,  ray  theory  is  not  valid  in  the  vicinity 
of  a  caustic  and  if  a  wave  theory  solution  were  utilized  in  this  re¬ 
gion  the  resultant  intensity  peak  would  probably  not  exceed  93  dB  or 
about  8  db  above  the  inverse-square-law  intensity  for  a  homogeneous 
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Fig.  — Intensity  variation  in  the  vicinity  of  the  caustic  ‘Case  4) 
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atmosphere.  Similarly  the  step  change  in  intensity  at  the  longer 

range  end  of  the  multiple  arrival  region  would  also  be  smoothed  out 

by  a  wave  solution  in  this  region.  However,  this  type  of  atmospheric 

profile  does  result  in  appreciable  deviations  from  an  inverse-square- 

law  drop-off  of  intensity  even  in  the  absence  of  atmospheric  attenua¬ 
tion. 


Case  5 

This  case  is  also  a  two-layer  atmosphere  with  characteristics  as 
specified  below. 
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If  as  in  Case  4  the  initial  wave  normal  azimuth  is  east  then 
Fig.  17  represents  the  plot  of  v  +  Wj(  for  this  azimuth  of  propagation 
A  comparison  of  this  figure  with  the  corresponding  one  for  Case  4 
shows  that  up  to  1000  ft  the  gradients  are  identical.  However,  above 

2000  ft  the  gradient  of  v  +  Wj(  decreases  to  a  value  of  7.5  ft/sec  per 
1000  ft. 

The  resulting  ray  pattern  for  a  source  at  500-ft  altitude  is 
shown  in  Fig.  18,  and  as  would  be  expected  from  the  positive  gradients 
of  Fig.  17,  the  rays  are  concave  toward  the  ground.  However,  the  cur¬ 
vature  of  those  portions  of  the  rays  above  1000  ft  is  considerably 
less  than  that  below  1000  ft  due  to  decrease  in  the  gradient  of  v  '+  w 
at  the  higher  altitude.  As  a  result,  the  ground  range  of  the  rays  in¬ 
creases  monotonically  with  initial  inclination  angle  and  there  is  no 
region  of  multiple  arrival  as  in  Case  4. 

In  Fig.  19  the  corresponding  variation  of  intensity  with  range 
shows  an  inverse-square-law  variation  out  to  19,361  ft.  This  range 
as  in  Case  4  corresponds  to  the  ray  which  is  tangent  to  the  1000  ft 
altitude  boundary  between  the  two  layers  of  the  assumed  atmosphere. 
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Beyond  this  critical  range,  the  intensity  drops  abruptly  to  zero,  af¬ 
ter  which  it  increases  rapidly  to  a  maximum  of  about  74.5  dB  at  32,500 
ft.  It  then  decreases  again,  always  lying  below  the  homogeneous  atmo¬ 
sphere  curve.  This  type  of  Intensity  variation  can  be  explained  as 
follows.  As  long  as  the  initial  elevation  angle  of  the  ray  is  low 
enough  that  the  ray  remains  below  the  1000-ft  level  the  divergence  of 
the  rays  results  in  the  inverse-square- lav  intensity  variation  observed 
previously  for  a  constant  gradient  of  v  +  w  .  However,  when  the  ini¬ 
tial  elevation  angle  reaches  a  value  at  which  the  ray  begins  to  pene¬ 
trate  the  upper  layer  with  a  lower  gradient  of  v  +  w^,  there  is  a  dis¬ 
continuity  in  the  rate  of  change  of  ground  range  with  elevation  angle. 

As  a  result  the  value  of  dr/d©^  in  Eq.  (113)  becomes  infinite  and  a 
zero  intensity  results  immediately  beyond  the  critical  range  of  19,361 
ft.  As  the  elevation  angle  continues  to  increase  the  value  of  dr/d© 

o 

decreases  and  the  intensity  increases  rapidly.  However,  as  the  eleva¬ 
tion  angle  increases  farther  and  more  and  more  of  the  path  lies  in  the 
upper  layer,  the  intensity  variation  begins  to  approach  the  inverse- 
square  law,  corresponding  again  to  a  constant  gradient  of  v  +  wx>  It 
should  also  be  noted  that  in  actual  practice  these  apparent  disconti¬ 
nuities  would  be  somewhat  smoothed  out  by  diffraction  effects  and  also 
by  the  fact  that  the  change  in  gradient  at  the  1000-ft  level  would  prob¬ 
ably  not  be  discontinuous  as  assumed  here. 

In  Fig.  20  the  intensity  variations  for  Cases  4  and  5  are  plotted 
on  the  same  graph  for  comparison.  It  is  seen  that  for  these  two  cases 
the  meteorological  conditions  up  to  1000  ft  altitude  are  identical. 
However,  as  a  result  of  differences  in  conditions  in  the  layer  between 
1000  and  2000  ft  there  are  differences  in  intensity  of  as  much  as  10 
dB  between  the  two  cases  at  ranges  in  the  vicinity  of  22,000  ft. 

Case  6 

The  assumed  meteorological  conditions  for  this  case  are  shown  be¬ 
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If  again  the  initial  wave  normal  azimuth  is  eastward,  then  Fig.  21 
represents  the  plot  of  v  +  Wj{  for  this  azimuth  of  propagation.  It  is 
seen  that  from  0  to  1000  ft  there  is  a  negative  gradient  of  v  +  w  of 
12.5  ft/sec  per  1000  ft,  while  from  1000  to  2000  ft  the  gradient  re¬ 
verses  and  has  a  value  of  37.5  ft/sec  per  1000  ft,  and  finally  above 

2000  ft  the  gradient  is  again  negative  with  a  value  of  12.5  ft/sec  per 
1000  ft. 


Figure  22  shows  the  resulting  ray  pattern  emanating  from  the 
source.  It  is  seen  that  in  the  regions  of  negative  gradient  of  v  +  w 
the  rays  are  concave  upward  while  in  regions  of  positive  gradient  they 
are  concave  downward.  For  initial  values  of  wave  normal  elevation  less 
than  -6.1  deg  the  value  (  f  C  is  greater  than  1100  ft/sec,  the  value  of 
v  +  wx  at  ground  level,  and  the  corresponding  rays  all  reach  the  ground 
The  ray  having  an  initial  wave  normal  elevation  of  -6.1  deg  has  a  value 
of  C  equal  to  1100  ft/sec  and  is  tangent  to  the  ground  at  a  range  of 
9192  ft,  after  which  it  rises  to  a  maximum  altitude  of  1333  ft,  where 
the  value  ofv+Wj(  again  equals  1100  ft/sec.  For  elevation  angles  be¬ 
tween  -6.1  deg  and  +6.1  deg  the  values  of  C  are  less  than  1100  ft/sec 
and  a  given  ray  is  bounded  between  the  two  altitudes  at  which  v  +  w 
equals  the  value  of  C  for  that  ray.  When  the  elevation  angle  equals 
+6.1  deg  the  ray  is  tangent  to  the  ground  at  a  range  of  25,871  ft  after 
first  reaching  a  maximum  altitude  of  1333  ft.  As  the  initial  elevation 
angle  increases  from  6.1  to  13.51  deg  the  values  of  C  lie  between  1100 
and  1125  ft/sec.  The  corresponding  rays  rise  to  a  maximum  altitude  at 
V  +  wx  e9ual  C  after  which  they  descend  to  ground  level.  As  the 
elevation  angle  increases  from  6.1  to  9.7  deg  the  ground  range  decrease; 
from  25,871  to  20,415  ft  while  from  9.7  to  13.51  deg  the  range  increase; 
again  to  21,617  ft.  For  elevation  angles  in  excess  of  13.51  deg  the 

value  of  C  is  greater  than  any  of  the  specified  values  of  v  +  w  and  th, 

x 


Alt*tude  (ft) 
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rays  do  not  reach  a  maximum  altitude  and  thus  do  not  return  to  ground 
level. 

This  example  demonstrates  the  phenomenon  of  ducting  since  all  of 
the  energy  radiated  between  elevation  angles  of  -6.1  deg  and  +6.1  deg 
is  confined  between  the  ground  and  an  altitude  of  1133  ft.  The  actual 
upper  bound  of  the  duct  is  the  nearly  horizontal  caustic  curve  shown 
in  the  figure.  A  second  caustic  curve  is  formed  by  the  rays  having 
initial  elevation  angles  between  6.21  and  9.7  deg.  This  caustic  is 
also  very  nearly  a  straight  line  and  intersects  the  ground  at  a  range 
of  20,415  ft.  There  are  two  shadow  zones  associated  with  this  ray 
pattern.  The  first  is  bounded  by  the  rising  portion  of  the  ray  tan¬ 
gent  to  the  ground  at  9192  ft,  the  second  caustic  curve,  ar.d  the 
ground  surface.  The  other  shadow  zone  is  entirely  above  ground  level 
and  is  bounded  by  the  two  branches  of  the  13.51  deg  ray  which  splits 
at  3000  ft  and  the  horizontal  caustic  which  is  the  upper  bound  of  the 
duct. 

In  Fig.  23  the  variation  of  intensity  with  ground  range  is  shown. 
It  is  seen  that  out  to  a  range  of  9192  ft  the  intensity  drops  off  in 
accordance  with  the  inverse-square  law.  Between  9192  and  20,415  ft 
the  intensity  drops  to  zero  while  from  20,415  to  25,871  ft  there  is 
a  second  zone  of  audibility  with  an  infinite  intensity  at  20,415  ft 
where  the  caustic  curve  intersects  the  ground,  and  a  region  of  double 
arrival  between  20,415  and  21,617  ft.  As  in  Case  4  diffraction  ef¬ 
fects  would  round  off  the  infinite  intensity  cusp;  however,  it  appears 
likely  that  in  this  second  region  of  audibility  the  intensity  might 
be  as  much  as  6  dB  above  the  curve  for  a  homogeneous  atmosphere. 

METEOROLOGICAL  DATA 

To  Investigate  the  nature  of  sound  transmission  under  arctic 
weather  conditions,  three  representative  sets  of  meteorological  data 
were  selected  from  Ref.  8.  The  information  relevant  to  sound  propa¬ 
gation  is  shown  in  Tables  1  through  3. 
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Table  1 
PROFILE  1 


Level 

Altitude 

(ft) 

Sound 

Ve loc ity 
(ft/sec) 

WE 

(ft/sec) 

WN 

(ft/sec) 

Relative 

Humidity 

(X) 

0 

0 

1010.8 

15.7 

8.9 

30.0 

1 

492 

1012.9 

30.2 

13.8 

35.0 

2 

820 

1015.0 

37.1 

13.8 

40.0 

3 

984 

1019.3 

40.7 

13.8 

43.0 

4 

2133 

1023.5 

52.2 

6.6 

57.0 

5 

4265 

1027. 7 

51.5 

-0.7 

53.0 

6 

6890 

1025.6 

49.2 

0.0 

41.0 

7 

9843 

1020. 3 

49.2 

0.0 

28.0 

Table  2 


PROFILE  2 


Level 

Altitude 

(ft) 

Sound 

Velocity 

(ft/sec) 

w 

E 

(ft/sec) 

W 

N 

(ft/sec) 

Relative 

Humidity 

(X) 

0 

0 

1050.7 

3.0 

2.6 

93.0 

1 

328 

1055.4 

19.4 

10.2 

93.0 

2 

656 

1060.3 

29.5 

8.9 

93.0 

3 

1148 

1067.0 

34.4 

3.3 

93.0 

4 

2625 

1069.0 

32.8 

0.0 

93.0 

5 

9843 

1044.5 

32.8 

0.0 

90.0 

Table  3 


PROFILE  3 


Level 

Altitude 

(ft) 

Sound 

Velocity 

(ft/sec) 

WE 

(f t/sec) 

WN 

(ft/sec) 

Relative 

Humidity 

(7.) 

0 

0 

1099.0 

36.7 

1 

79.0 

1 

656 

1096.2 

47.6 

79.5 

2 

1640 

1092.1 

63.0 

80.0 

3 

2953 

1086.5 

68.6 

81.0 

4 

4265 

1081.1 

68.2 

0.0 

82.5 

5 

7874 

1067.0 

66.6 

0.0 

85.0 

6 

9843 

1061.9 

65.6 

0.0 

85.0 
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ARCTIC  SOUND  TRANSMT^mM 

In  the  following  paragraphs,  the  variation  of  sound  intensity  as 
a  function  of  ground  range  and  azimuth  of  propagation  is  determined 
for  each  of  the  meteorological  profiles  given  in  Tables  1  through  3. 

In  all  of  these  cases  a  source  power  of  100  kW  at  an  altitude  of  500 
ft  is  assumed. 

Profile  1 

In  Fig.  24  plots  of  v  +  as  a  function  of  altitude  are  shown 
for  Profile  1  at  intervals  of  30  deg  in  azimuth  from  0  to  360  deg, 
where  the  azimuth  angle  is  measured  counterclockwise  from  east. 

The  variation  of  intensity  as  a  function  of  ground  range  is  shown 
in  Fig.  25  for  an  azimuth  angle  of  0  deg.  It  is  seen  that  out  to  a 
range  of  11,500  ft  the  intensity  has  essentially  an  inverse-square- law 
decrease.  Between  11,500  and  13,000  ft  an  increase  occurs  because  of 
the  increase  in  the  gradient  of  v  +  w^  at  the  820-ft  level.  This  ef¬ 
fect  is  similar  to  that  shown  in  Example  4  and  results  in  a  caustic 
which  intersects  the  ground  at  about  11,750  ft  range.  The  branch  of 
the  intensity  curve  between  12,000  and  34,000  ft  range  is  due  to  the 
decrease  in  gradient  at  the  984-ft  level  and  is  similar  to  Example  5. 
Finally,  the  branch  beyond  34,000  ft  in  range  results  from  the  decrease 
in  gradient  at  the  2133-ft  level. 

As  the  azimuth  angle  increases,  the  gradient  in  the  layer  above 
2133  ft  becomes  negative  and  the  corresponding  branch  of  the  intensity 
curve  disappears  as  shown  in  Fig.  26  since  rays  whi.h  go  above  this 
level  are  no  longer  refracted  back  to  the  ground.  A  further  increase 
in  azimuth  results  in  a  negative  gradient  above  980  ft  so  that  the 
second  branch  of  the  intensity  curve  also  disappears  as  shown  in  Fig. 

27.  For  azimuth  angles  between  120  and  210  deg  the  gradient  is  com¬ 
pletely  negative  and  the  intensity  curve  is  similar  to  that  shown  for 
the  upwird  direction  in  Example  3  where  a  shadow  zone  is  formed. 

At  an  azimuth  of  240  deg  the  shadow  zone  still  exists,  however, 
because  of  the  positive  gradient  above  820  ft  the  value  of  v  +  w 
above  3000  ft  exceeds  that  at  ground  level,  and  it  is  possible  to  have 
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rays  which  rise  initially  to  altitudes  between  3000  and  4000  ft  before 

returning  to  ground  level  at  ranges  of  the  order  of  70,000  ft.  The 

resulting  intensity  variation  is  shewn  in  Fig.  28.  Finally,  for  azi¬ 
muths  beyond  about  270  deg  the  intensity  curves  again  take  on  the  form 
shown  in  Fig.  25. 

Profile  2 

In  Fig.  29  the  plots  of  v  +  w^  at  30-deg  increments  in  azimuth 

are  shown  for  Profile  2.  For  azimuths  between  0  and  90  deg  the  inten¬ 

sity  variation  is  similar  to  that  shown  in  Example  5  and  results  from 
the  decrease  in  gradient  at  the  656-ft  level.  Figure  30  is  typical  of 
this  range  of  azimuth  angles. 

For  azimuth  angles  from  about  120  deg  up  to  210  deg  the  resulting 
intensity  curves  are  again  similar  to  Example  3  in  the  upwind  direc¬ 
tion,  with  a  shadow  zone  at  ranges  beyond  5000  to  8000  ft. 

For  azimuths  between  240  and  270  deg  the  higher  value  of  v  +  w 
at  altitudes  between  1000  and  2500  ft  results  in  rays  being  refracted 
back  to  the  ground  from  these  altitudes  at  ranges  in  excess  of  20,000 
ft. 

Finally,  for  azimuths  between  300  and  360  deg  the  intensity  curves 
are  similar  to  that  shown  in  Fig.  30  but  with  an  additional  low  inten¬ 
sity  branch  beyond  16,500  ft  range  resulting  from  transmission  through 
the  very  slight  positive  gradient  between  1148  and  2625  ft,  as  shown  in 
Fig.  31. 

Profile  3 

For  Profile  3  the  altitude  dependence  of  v  +  w^  at  30-deg  incre¬ 
ments  in  azimuth  ia  shown  in  Fig,  32. 

For  azimuth  angles  between  0  and  60  deg  the  intensity  variation 
is  similar  to  that  shown  in  Example  5  due  to  the  decrease  in  gradient 
at  the  656-ft  level.  Figure  33  shows  the  variation  for  60-deg  azimuth. 
The  maximum  range  of  41,289  ft  corresponds  to  a  ray  with  a  maximum  al¬ 
titude  of  1640  ft,  and  any  ray  which  exceeds  this  altitude  is  in  a  re¬ 
gion  of  negative  gradient  and  cannot  return  to  the  ground  level. 
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32  —  Dependence  of  v  +  on  oltitude  at  30° 
in  azimuth  of  propagation  (Profile  3) 
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For  azimuths  between  90  and  270  deg  the  gradient  is  negative  at 
all  altitudes  and  a  shadow  zone  is  formed  similar  to  that  in  Example  3. 

Finally,  for  azimuths  between  270  and  360  deg  there  is  a  possibil¬ 
ity  of  an  enhancement  of  intensity  similar  to  that  in  Example  4  due  to 
the  increase  in  gradient  at  the  656-ft  level  and  a  reduction  of  inten¬ 
sity  similar  to  that  in  Example  5  due  to  the  decrease  in  gradient  at 
1640  ft.  This  is  illustrated  in  Fig.  34. 


ATTENUATION  EFFECTS 

It  is  found  that  the  classical  attenuation  due  to  viscosity  and 
heat  conduction  is  negligible  in  comparison  with  that  due  to  molecular 
absorption.  In  Tables  4,  5,  and  6  the  average  temperature  and  humidity 
m  each  atmospheric  layer  and  the  molecular  absorption  coefficients  for 
frequencies  of  125,  250,  500,  and  1000  cps  are  presented  for  each  of 
the  meteorological  profiles. 

An  examination  of  Table  4  shows  that  for  the  low  water  content  of 
the  atmosphere  at  these  temperatures,  the  resulting  molecular  attenua¬ 
tion  is  less  than  0.1  dB  per  1000  ft  even  at  a  frequency  of  1000  cps. 

Thus  the  variation  of  intensity  with  distance  is  essentially  indepen¬ 
dent  of  frequency  for  Profile  1. 

A  similar  examination  of  Table  6  shows  that  the  attenuation  in¬ 
creases  with  both  frequency  and  altitude.  However,  since  the  sound 
rays  which  reach  the  ground  from  the  assumed  500-ft  source  altitude 
never  reach  layers  4,  5,  or  6,  the  higher  values  of  attenuation  are  not 
encountered.  Figure  35  shows  the  variation  of  intensity  with  distance 
for  frequencies  of  125  and  1000  cps  for  the  attenuation  shown  in  Table 
6.  In  this  case  the  sound  rays  are  confined  to  the  first  two  layers  of 
the  atmosphere  where  the  attenuation  of  the  125-cps  frequency  is  about 
0*01  dB  Per  1000  ft  while  that  for  1000  cps  is  from  0.35  to  0.48  dB  per 
1000  ft.  Thus  at  the  longer  ranges  the  1000-cps  curve  lies  10  to  15  dB 
below  that  for  the  lower  frequency. 

In  the  case  of  Profile  2  as  shown  in  Table  5,  the  combination  of 
temperature  and  humidity  results  in  comparatively  high  attenuation  at 
all  altitudes,  ranging  from  about  0.15  dB/1000  ft  at  125  cps  to  about 
4  dB/1000  ft  at  1000  cps.  Figure  36  shows  an  example  of  the  variation 
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Fig.  34 — Intensity  variation  with  range  at  330°  azimuth  (Profile  3) 


-73- 


Table  4 


ATMOSPHERIC  ATTENUATION  FOR  PROFILE  1 


Relative 

Temper- 

Attenua  t ion 

Coefficient 

Layer 

Humid ity 

ature 

(dB/kft) 

(7.) 

(deg  C) 

125  cps 

250  cps 

500  cps 

1000  cps 

1 

32.5 

-36.6 

0.0077 

0.0077 

0.0077 

0.0077 

2 

37.5 

-35.6 

0.0129 

0.0129 

0.0129 

0.0129 

3 

41.5 

-34.  1 

0.0221 

0.0221 

0.0221 

0.0221 

4 

50.0 

-32. 1 

0.0499 

0.0502 

0.0502 

0.0502 

5 

55.0 

-30. 1 

0.0933 

0.0946 

0.0949 

0.0950 

6 

47.0 

-29.6 

0.0767 

0.0774 

0.0776 

0.0776 

7, 

34.5 

-31.4 

0.0282 

0.0283 

0.0283 

0.0283 

Table  5 

ATMOSPHERIC  ATTENUATION  FOR  PROFILE  2 


Relative 

Temper- 

Attenuation 

Coefficient 

Layer 

Humidity 

ature 

(dB/kft) 

(4) 

(deg  C) 

125  cps 

250  cps 

500  cps 

1000  cps 

1 

93.0 

-17.0 

0. 1831 

0.6601 

1.8927 

3.5498 

2 

93.0 

-14.6 

0.1250 

0.4787 

1.6353 

4.1298 

3 

93.0 

-11.8 

0.0773 

0.3045 

1.1511 

3.7739 

4 

93.0 

-9.6 

0.0534 

0.2122 

0.8278 

3.0135 

5 

91.5 

-15.2 

0.1409 

0.5325 

1.7447 

4.0492 

Table  6 

ATMOSPHERIC  ATTENUATION  FOR  PROFILE  3 


Layer 

Relative 

Humidity 

(7.) 

— - - 1 

Temper¬ 
ature 
(deg  C) 

Attenuation  Coefficient 
(dB/kft) 

125  cps 

250  cps 

500  cps 

1000  cps 

1 

79.3 

5.1 

0.0056 

0.0222 

0.0889 

0.3552 

2 

79.8 

3.4 

0.0075 

0.0298 

0.1194 

0.4769 

3 

30.5 

0.9 

0.0113 

0.0450 

0. 1800 

0.7180 

4 

81.8 

-1.8 

0.0177 

0.0706 

0.2820 

1.1191 

5 

83.8 

-6.7 

0.0392 

0.1565 

0.6185 

2.3619 

6 

85.0 

-11.4 

0.0862 

0.3387 

1.2649 

3.9991 
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of  intensity  with  distance  for  frequencies  of  125,  250,  500  and  1000  cp 
It  is  seen  that  even  at  a  range  of  8500  ft  there  is  a  difference  of 
30  dB  between  the  125-cps  curve  and  that  for  1000  cps. 


IV.  CONCLUS IONS 


On  the  basis  of  the  foregoing  results  the  following  specific  con¬ 
clusions  can  be  stated  regarding  sound  propagation  from  an  elevated 
source  to  the  ground. 

o  A  sound  ray  emanating  from  an  elevated  source  with  an  initial 
inclination  angle  less  than  or  equal  to  zero  will  reach  ground 
level  if  its  refraction  constant  C  (see  Eq.  (11))  is  greater 
than  the  maximum  value  of  v  +  Wj(  existing  between  the  source 
and  the  ground. 

o  A  sound  ray  with  an  initial  inclination  greater  than  zero  will 
reach  ground  level  if  its  value  of  C  not  only  satisfies  the 
previous  condition  but  also  is  less  than  the  maximum  value  ol' 
v  +  wx  occurring  above  the  source. 

o  If  a  ray  has  an  initial  elevation  angle  less  than  zero  and  its 

value  of  C  is  equal  to  the  maximum  value  of  v  +  w  existing  fc#-- 

x 

low  the  source,  then  its  ground  range  represents  the  boundary 
of  a  shadow  zone.  Propagation  of  acoustic  energy  beyond  this 
boundary  can  occur  only  by  scattering  and  diffraction  and  its 
intensity  drops  off  quite  rapidly  with  distance  beyond  the 
boundary. 

o  If  the  maximum  value  of  v  +  w^  as  a  function  of  altitude  occurs 
above  the  source  then  a  ray  with  a  positive  elevation  angle  and 
a  value  of  C  equal  to  this  maximum  has  a  ground  range  which 
represents  the  maximum  range  of  direct  acoustical  propagation. 

o  Variations  of  v  +  w^  between  source  altitude  and  the  ground 
appear  to  have  relatively  little  effect  cn  the  intensity  at 
ground  level.  Thus  in  that  area  of  the  ground  illuminated  by 
rays  whose  initial  elevation  angles  are  negative  the  intensity 

variation  with  distance  deviates  very  little  from  the  inverse- 
square  law. 

o  If  the  gradient  of  v  +  w^  decreases  d iscontinuously  at  some 
altitude  above  the  source,  then  the  intensity  at  ground  level 
in  the  vicinity  of  those  rays  which  just  penetrate  the  upper 
layer  will  be  abnormally  low. 
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o  If  the  gradient  of  v  +  wx  increases  discontinuously  at  some 
altitude  above  the  source,  the  intensity  at  ground  level  in 
the  vicinity  of  those  rays  which  just  penetrate  the  upper 
layer  will  be  abnormally  high. 

o  Atmospheric  attenuation  reduces  the  intensity  of  the  higher 
frequency  components  significantly  but  has  comparatively  lit¬ 
tle  effect  at  frequencies  of  the  order  of  125  cps,  which  are 
readily  available  in  the  spectrum  of  multi-engine  jet  aircraft, 
o  For  the  assumed  source  altitude  of  500  ft  and  for  the  three 
meteorological  conditions  selected  the  source  could  be  heard 
at  ground  level  out  to  at  least  6000  ft  in  any  direction.  How¬ 
ever,  this  value  would  increase  with  source  altitude  and  with 
source  power  and  would  vary  with  azimuth  of  propagation  depend¬ 
ing  on  tie  wind  structure.  As  a  result  It  is  extremely  diffi¬ 
cult  to  specify  a  detection  range  for  a  ground-based  acoustic 
sensor  used  to  detect  elevated  targets. 
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